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1.  INTRODUCTION 


This  annual  report  covers  the  second  year  of  fontract  Nc. 

F 19628-8 3 -K- 9028 .  It  consists  of  two  submitted  :apers,  one  to  a 
Conference  and  one  to  a  Journal. 

Section  2  of  this  report  is  a  paper  presented  at  the  NATO  AJAR 
Spring  1986  Electromagnetic  Wave  Propagation  Panel  Symposium  on  "Tl 
Aerospace  Environment  at  High  Altitudes  and  its  Imp lications  for 
Spacecraft  Charging  and  Communications",  The  Hague,  Netherlands, 
2-6  June  1986.  This  paper  is  primarily  of  a  "survey"  nature,  and 
contains  relatively  brief  accounts  of  several  aspects  of  the 
low-p.olar-orbit  high-voltage  charging  problem.  Howe'er,  its  Secti 
presents  the  first  results  from  L.w.  Parker’s  simulation  work  whig' 
was  supported  by  this  Contract.  We  believe  that  his  study  is  the 
first  calculation  of  downstream  potentials  on  electrically-isolate 
surfaces  which  is  based  on  the  use  of  ion  orbit  integrations  for 
calculating  current  deposition  on  such  surfaces.  The  successful 
completion  of  this  work  involved  overcoming  extreme  numerical  dir: 
These  appear  to  be  closely  associated  with  physical  sensitivities 
inherent  in  the  physics  of  high-voltage  ion  wakes,  particularly  i 
regard  to  the  extremely  beam-like  character  of  the  ion  'elocity 
distributions  in  such  wakes.  This  situation  has  severe  implicatio 
for  the  ability  to  make  economical,  realistic  predictions  or 
high-voltage  charging  in  geometrical ly-conpl icated  situations, 
including  arrangements  of  exposed  equipment  boxes  m  the  Shuttle's 
cargo  bay,  and  much  of  the  Space  Station's  proposed  Polar  Platrcrm 


Section  3  of  this  report  contains  the  text  of  a  paper  submitted 
in  October  1936  for  publication  in  Journal  of  Geophysical  Research, 
Space  Physics.  It  is  an  expanded  version  of  an  earlier  paper  by 
baf ranbcise  (1935)  on  secondary-electron  escape  from  negatively-charge- 
spacecraft  surfaces  in  a  magnetic  field,  and  now  includes  results  for 
situations  in  which  there  exist  electric  fields  tangential  as  well 
as  normal  to  such  surfaces. 
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SPACECRAFT  CHARGING  IN  THE  AURORAL  PLASMA: 
PROGRESS  TOWARD  UNDERSTANDING  THE  PHYSICAL  EFFECTS  INVOLVED 


J.G.  Laframboise*  and  L.W.  Parker** 

•Physics  Department,  York  University,  Toronto,  Canada  M3J  1P3 
••Lee  W.  Parker,  Inc.,  252  Lexington  Rd.,  Concord,  Mass.,  U.S.A.,  017-12 

ABSTRACT 


The  work  presented  here  is  in  four  parts.  In  the  first,  we  review  the  main 
differences  between  the  oiasma  environments  in  geostationary  orbit  and  low  poiar 
orbit  with  regard  to  high-voltage  charging  situations.  We  next  present  results  from  a 
calculation  of  secondary-electron  escape  currents  from  negatively-charged  spacecraft 
surfaces  having  various  orientations  relative  to  the  local  magnetic-field  direction. 

We  show  that  for  finite  ranges  of  combinations  of  electric  ana  magnetic  field  direc¬ 
tions,  secondary-electron  escape  is  completely  suppressed  and  therefore  cannot  help 
to  discharge  the  spacecraft.  In  such  circumstances,  secondary  electrons  may  travel 
distances  many  times  their  gyro  radii  before  reimpacting,  and  this  may  produce 
greatly  increased  secondary-electron  surface  currents.  Thirdly,  we  develop  a  simple 
rough  estimate  of  the  required  conditions  for  high-voltage  auroral-zone  charging,  ihe 
resuits  suggest  that  for  any  given  spacecraft,  surface  potentials  are  likely  to  depend 
more  strongly  on  the  ratio  or  ambient  flux  of  high-energy  electrons  to  that  of  all  ions, 
than  on  any  other  environmental  parameter.  Finally,  we  present  preliminary  results 
of  numerical  simulation  work  directed  toward  testing  this  hypothesis.  Numerical 
instabilities  enccimtered  in  doing  this  simulation  work  probably  are  closely  related  to 
physical  sensitivities  inherent  in  the  physics  of  the  ion  wake  behind  the  soacecraft, 
and  especially  to  beam-iike  constituents  of  the  ion  population  in  the  wake. 


1.  INTRODUCTION 

The  plasma  environment  in  low  Earth  orbit  has  very  different  properties  from  that 
in  geostationary  orbit.  In  GEOSTATIONARY  ORBIT: 

(a)  the  Debye  Length  is  a  few  tens  of  metres;  therefore,  SPACE-CHARGE 
COUPLING  is  SMAiJ_,  and  one  can  ignore  it  entirely,  or  use  a  linearized  approx¬ 
imation  for  it,  wnen  calculating  potentials  near  soacecraft  of  ’’ordinary''  size. 

(b)  average  particle  gyrorariii  are  a  few  tens  of  metres  or  larger;  therefore 
MAGNETIC-FIELD  Elrrc.uiS  ARc.  NEGLIGIBLE  for  calculating  particle  orbits,  near 
spacecraft  of  "ordinary"  size. 

!c)  circular-orbit  soeed  is  much  smaller  than  particle  average  random  speeds; 
therefore.  AMBIENT  r  ARTICLE  DISTRIBUTIONS  CAN  “USUALLY"  BE  ASSUMED 
ISOTROPIC,  and  this  greatly  simplifies  calculations  of  ambient-particle  currents 
collected  by  spacecraft.  In  reality,  significant  anisotropies  are  observed  in  ambient 
particie  distributions,  out  these  are  generally  less  important  than  the  relative  ion 
anisotropy  induced  in  low  orbit  by  spacecraft  motion. 

.n  LCW  EARTH  CR3IT,  none  of  the  above  are  true.  Space-charge  coupling  is  strong, 
magnetic-: ;e.d  effects  are  strong,  and  ion  flow  effects  are  important.  Also,  charged- 
particie  mean  Tree  paths  may  not  aiwavs  be  much  larger  than  spacecraft  size, 
especially  during  thruster  firings  and  water  dumps. 

.ac.e  1  summarizes  some  important  characteristic  lengths  and  soeeds  for  low- 
orbit  conditions.  A  surprising  feature  of  this  Table  is  that  the  sheath  thicknesses 
indicated  are  much  larger  than  the  ambient  Debye  Lengtn,  but  this  is  because  the 
sheath  potentials  are  much  larger  than  the  ambient-particle  thermal  energies.  These 
distances  are  at  most  crmparaoie  to  typical  spacecraft  dimensions,  in  contrast  with 
‘.he  geosynchronous  situation.  Table  1  also  shows  that  secondary  electrons  nave  an 
average  gyroraoius<<  typical  soacecraft  dimensions,  so  their  escaoe  will  be  inhibited 
strongly  on  surfaces  which  are  nearly  paranei  to  the  magnetic  field  B  (Fig.  1),  while 
aurorai  electrons  have  an  average  jjyroradius  £  typical  spacecraft  dimensions,  so 
their  collection  will  be  affected  omy  moderately,  exceot  for  very  large  spacecraft. 

>ve  return  to  this  Question  in  Sec.  2.  Also  evident  from  Table  1  is  the  large  value  of 
t*ve  ion  speed  ratio  (spacecraft  sceed/ion  most-probable  thermal  speed)  in  low-orbit 
concit.ons.  In  these  conditions,  ion  collection  on  oownstream  surfaces  will  be 
.mcited.  If  a  surface  is  simultaneously  oownstream  and  nearly  parallel  to  the 
—  apnetic  fie.c,  as  is  Lkeiy  to  be  the  case  in  the  aurora i  zones,  then  the  tendency  for 
nigh-voitage  charging  to  occur  on  it  will  be  greatly  increased  (Fig.  2). 
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In  the  rest  cf  this  paper,  we  report  briefly  on  three  separate  projects.  _  In  Seq.  2, 
we  present  resuits  of  a  calculation  of  secondary-electron  escape  currents  from 
negatively-charged  surfaces  having  various  orientations  relative  to  the  local  mag¬ 
netic-field  direction.  In  Sec.  3,  we  develop  a  simple  rough  estimate  of  the  required 
conditions  for  high-voltage  charging.  The  results  suggest  that  for  any  given  space¬ 
craft,  surface  potentials  are  likely  to  deoend  more  strongly  on  the  ratio  of  ambient 
flux  of  high-energy  electrons  to  that  of  all  ions,  than  on  any  other  environmental 
parameter.  In  Sec.  4,  we  present  results  of  numerical  simulation  work  directed 
toward  testing  this  hypothesis.  This  work  involves  calculations  of  floating-potential 
distributions  on  infinite  cylinders  in  collisionless  plasma  crossflows  whose  prop¬ 
erties  model  those  of  the  auroral  plasma. 

2.  CALCULATION  OF  SECONDARY-ELECTRON  ESCAPE  CURRENTS  FROM 
NEGATIVELY-CHARGED  SURFACES  IN  A  MAGNETIC  FIELD. 

Figure  2  illustrates  why  escape  of  secondary  electrons  Is  affected  by  magnetic 
fields.  In  Fig.  2(a),  the  spacecraft  surface  is  perpendicular  to  B,  and  the  emitted 
electrons,  which  experience  an  electric  force  -eE  directed  away  from  the  surface, 
all  escape,  helping  to  discharge  it.  Here  e  is  the  magnitude  of  the  electronic  charge. 
In  Fig.  2(b),  the  spacecraft  surface  is  nearly  parallel  to  B,  and  almost  all  of  the 
emitted  electrons  return  to  it,  even  though  they  still  experience  an  electric  force 
directed  away  from  it.  These  electrons  therefore  are  inable  to  help  discharge  it,  so 
a  surface  nearly  parallel  to  B  Is  more  likely  to  charge  to  a  large  negative  voltage. 

The  component  or  E  which  is  perpendicular  to  B  results  oniy  in  an  E  x  B  drift  parallel 
to  the  surface. 

In  this  Section,  we  present  numerically-calculated  escaping  secondary-electron 
fluxes  for  these  conditions,  for  the  case  when  the  eiectric  field  E  is  normal  to  the 
surface.  We  aiso  indicate  some  general  properties  of  the  escaping  flux  when  E  has  a 
nonzero  tangential  component,  l.e.,  the  surface  is  at  a  noniniform  potential. 

We  assume  that  the  surface  Is  flat,  that  E  and  B  are  uniform,  and  that  secondary 
electrons  are  emitted  with  a  Maxwellian  velocity  distribution  corresponding  to  a 
temperatur  ->  T.  If  -eE  is  directed  along  the  outward  surface  normal,  then  the  ratio 
i  =  ;/!  of  escaping  to  emitted  flux  is  a  function  of  only  two  parameters:  the  angle 

0  between  the  surface  normal  and  the  direction  of  B  (Fig.  3),  and  a  parameter 
describing  the  strength  of  E.  A  convenient  choice  for  this  parameter  is  the 
difference  .in  potential  across  a  mean  secondary-electron  gyro  radius  a  =  (1/eB) 
(n-rr.kT/2)?,  divided  by  kT/e,  where  m  is  electron  mass  and  k  Is  Boltzmann’s 
constant.  This  quotient  is: 

,  _  E  /  irm 

where  E  =  IE|  and  S  s  |B(.  (2.1) 

This  quantity  aiso  has  an  alternative,  more  useful  interpretation:  it  is  the  ratio 
of  the  magnitude.  (E  x  Bi/B2  of  the  E  x  B  drift  speed,  to  one-half  the  mean  thermal 
speed  (SkT/irm)*  of  the  emitted  electrons. 

Our  method  of  calculating  escaping  flux  is  as  follows.  For  each  of  a  sufficient 
number  of  values  of  0  and  e,  we  cnoose  a  large  enough  discrete  set  of  values  of 
emission  velocity  components  vxc,  v  ,  and  vzq;  we  then  search  numerically  along 

each  of  the  resulting  electron  orbits,  which  are  known  analytically  with  time  as  a 
parameter,  in  order  to  determine  which  ones  reimpact  the  surface.  We  then  sum 
over  the  escaping  ones  to  calculate  the  flux,  and  we  tabulate  the  resulting  flux  values 
as  a  function  of  8  and  «.  Further  details  of  the  calculation  have  been  given  by 
Laframboise  (1935).  The  resulting  values  of  escaping  secondary-electron  current 
density  a re  shown  in  Fig.  4.  These  results  are  accurate  to  within  about  0.5%  or 
better  (Laframboise,  1985’.  An  empir'cal  analytic  expression  which  aoprcximates 
these  resuits  to  within  3%  of  J  has  been  given  by  Laframboise  11985,  £q.  3.1). 

An  imoortar.t  feature  of  these  results  is  that  when  f  is  large  enough,  electron 
escaoe  becomes  essentially  complete  except  when  0  is  very  nearly  90°.  This  means 
that  m  Shuttle  r.igh-voitage  charging  conditions,  for  which  30  S  (  i  120 
Laframboise,  1965),  the  occurrence  of  high-voltage  charging  in  marginal 
circumstances  may  depend  very  strongly  on  the  precise  orientation  cf  a  surface.  A 
sicwiv-rotatmg  surface  which  Dasses  through  tangentiahty  to  B  may  experience  a 
sudden,  brief  high-voltage  charging  event.  For  the  same  reason,  attempts  to  predict 
hignvoitage  charging  mav  be  afflicted  by  sensitivity'  problems:  if  one  attempts  to 
precict  worst -case  charging  by  assuming  that  secondaries  do  not  escape,  then  the 
resulting  predictions  are  likely  to  be  overly  pessimistic  most  of  the  time.  On  the 
utner  hand,  if  one  assumes  that  secondaries  do  escape,  correct  predictions  will  be 
obtained  almost  ali  cf  the  time,  but  occasionally  a  large  underestimate  of  charging 
will  occur. 


m 
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If  a  spacecraft  surface  Is  charged  to  a  nonuniform  potential,  then  E  will  have  a 
nonzero  component  tangential  to  the  surface.  The  resulting  problem  geometry  Is 
more  complicated,  and  two  additional  angles,  a  and  <//,  whicn  define  tne  direction  of 
-eE,  need  to  be  specified  (Fig.  5).  Calculations  of  escaping  fluxes,  and  also  the 
surface  currents  produced  by  retmpactlng  electrons,  are  in  progress  (J.G. 
Laframboise,  to  be  published).  Here  we  indicate  only  some  general  features  which 
these  results  must  possess. 

It  may  happen  that  even  though  -e£  points  outward  from  the  surface.  Its  projection 
along  B  points  inward  (Fig.  6a).  This  will  happen  If  the  angle  between  -eE  and  B  (or 
-B,  if  this  points  outward)  is  larger  than  90°.  In  this  case,  every  electron 
experiences  an  average  acceleration  toward  the  surface,  and  no  electrons  escape. 

This  will  greatly  enlarge  the  range  of  surface  orientations  for  which  electron  escape 
is  strongly  inhibited;  wnen  -eE  was  normal  to  the  surface.  Fig.  4  indicated  that  this 
range  was  a  few  degress  or  less.  The  situation  Is  analogous  to  that  for  a  sailboat 
tacking  against  the  wind  (Fig.  6b). 

The  situation  in  which  all  electrons  reimpact  can  be  further  subdivided  (Fig.  7). 
The  z  -component  of  the  E  x  B  drift  direction  is  -BxE ^/B2.  If  Bx  >  0 ,  this  is  opposite 

in  sign  to  E^.  Therefore,  if  E^  <  0 ,  electrons  which  do  not  reimpact  during  the 

first  gyropenod  after  their  emission  are  likely  to  travel  much  farther  along  the 
surface  before  they  do  reimpact,  and  the  resulting  surface  current  is  likely  to  be 
much  larger,  than  if  E^  >  0  (Fig.  7).  Another  consequence  If  E^  <  0  is  that 

electrons  are  mere  likeiy  to  travel  distances  across  which  our  assumption  of 
uniformity  of  E  and  B  is  no  longer  valid,  and  our  detailed  predictions  will  then 
become  inaccurate. 

3.  ESTIMATE  CF  REQUIRED  ENVIRONMENTAL  CONDITIONS  FOR  LOW-POLAR- 
CRBIT  CHARGING. 

In  this  Section,  we  develop  a  simple  rough  estimate  of  the  conditions  necessary 
for  high-voitage  charging  in  iow  polar  orbit,  and  we  show  that  spacecraft  surface 
potentials  are  iikely  to  depend  more  strongly  on  the  ratio  of  ambient  flux  of  high- 
energy  electrons  to  that  of  al^  ions  than  on  any  other  applicable  environmental 
parameter.  To  do  this,  we  make  the  following  approximations: 

(!)  We  assume  that  magnetic- field  effects  on  charged-particle  motion  are  negli¬ 
gible.  This  assumption  should  be  acceptable  for  initial  estimates  because  the  gyro- 
ratm  of  ions  and  nigh-energy  electrons  are  generally  a  few  metres  or  larger,  espe¬ 
cially  in  a  high-voitage  sheath  (Table  1),  and  collection  of  "cold"  ('0.1  eV)  iono¬ 
spheric  electrons  bv  a  negativeiy-charged  spacecraft  will  be  very  small,  so  their 
density  is  weil-approximated  by  a  Boltzmann  factor,  independently  of  the  presence  of 
a  magnetic  field. 

(2)  We  assume  that  ambient  high-energy  electrons  have  an  Isotropic  velocity 
dlstrioution.  Large  departures  from  this  have  been  observed  in  auroral -plasma 
conditions  (W.J.  Burke,  1984,  private  communication),  but  this  should  not  seriously 
affect  the  type  of  rough  estimate  made  here.  Parks  and  Katz,  1981,  and  Katz  and 
Parks,  1983,  assumed  both  the  ion  and  electron  fluxes  to  be  unidirectional;  we 
discuss  this  point  later  in  this  Section. 

(3)  We  ignore  secondary-electron  emission;  magnetic-field  effects  would  tend  to 
suppress  irus  on  some  parts  of  the  spacecraft  in  any  case  (Laframboise,  1983a, 

1985;  Sec.  2). 

(4)  We  assume  that  the  spacecraft  is  a  unpotentlal  sphere,  large  compared  to  the 
typical  ambient  Debve  length  of  £  1  cm.  We  consider  only  overall  charging  of  the 
spacecraft.  :his  regiects  the  possibility  that  local  high-voitage  charging  may  occur, 
especially  on  surfaces  in  the  spacecraft  wake  (Sec.  4). 

15)  We  assume  that  both  ions  and  electrons  have  double-Maxwell lan  velocity  distri¬ 
butions,  with  the  colder  component  in  either  case  having  a  temperature  of  0.1  eV, 
ana  the  hotter  1  keV  or  larger .  In  the  spacecraft  reference  frame,  these  are  super¬ 
posed  on  a  drift  velocity  equal  and  opposite  to  the  spacecraft  velocity. 

(6)  Ions  are  assumed  to  be  either  H*  or  O*  . 

Note  that  assumption  ;3)  could  cause  a  false  prediction  that  high-voltage  charging 
occurs,  while  assumotion  i4)^couid  cause  a  false  prediction  that  it  does  not.  The 
effects  of  assumptions  (1),  (2),  and  (5)  are  less  clear;  these  could  conceivably  either 
Increase  or  decrease  predicted  surface  potentials.  With  regard  to  (6),  assuming  that 
the  ions  are  H*  results  in  maximum  wake-filling  by  ions,  .t  there  are  any  elec'ri- 
caiiy-isoiatea  surfaces  in  the  spacecraft  wake,  tms  would  result  m  decreased  surface 
potentials  (magnitudes);  assuming  O’  gives  the  reverse. 


r»!i 


Probably  the  most  serious  difficulty  In  formulating  a  theory  for  low-orbit 
charging  Is  the  prediction  of  Ion  collection  on  downstream  surfaces.  As  mentioned  in 
assumption  (3)  above,  we  avoid  this  difficulty  by  considering  only  total,  rather  than 
local,  ton  collection,  on  a  unipotenttal  sphere.  Kanal  [1962,  Eq.  (63)]  gives  an 
expression  for  the  ion  current  collected  by  such  a  sphere  from  a  drifting  Maxwellian 
plasma  in  the  limit  of  zero  potentials  (relative  to  space  potential),  as  fallows: 

lj  =  i  l*1  (S1  +  -25-)  erftSj)  +  expt-Sj*)]  (3.1) 

where  1,  =  lj/1^,  IQl  Is  the  ion  random  current  en^kT^irm^,  S,  = 

U/(2kTj/mj)^  is  the  ion  speed  ratio,  U  is  the  ion  drift  speed  relative  to  the  space¬ 
craft,  e  is  the  magnitude  of  the  electronic  charge,  k  is  Boltzmann’s  constant,  and  nr, 
T.,  and  niaj  are  ion  mass,  temperature,  and  ambient  number  density.  We  assume 
that  U  =  8  km/sec,  corresponding  to  low  circular  orbit. 

We  need  to  take  account  of  the  efFect  of  a  large  ion-attracting  surface  potential  on 
ion  collection,  in  the  limit  of  small  Debye  length  Xp  compared  to  the  sphere  radius 

rg.  To  do  this,  we  use  a  result  of  Parrot  et  al  (1982).  These  authors  show  that  for 
a  probe  in  a  collisionless,  nonmagnetized,  Maxwellian  plasma  having  T^/T  =  1  and 
without  ion  drift,  and  in  the  limit  when  XQ/rs  — ►  0  but  -ed>s/kT  >>  1  [where  $s  is 

surface  potential  relative  to  space,  and  these  limits  must  be  approached  in  such  a 
way  that  (-e^s/kT) (Xj-j/rs) 4/3  remains  <<  1,  i.e.,  sheath  thickness  remains  <<  sphere 

radius],  the  ion  (attracted-particle)  current  is  larger  than  the  random  current  by  a 
factor  of  1.45.  This  factor  represents  the  effect  of  "presheath"  electric  fields  on  ion 
collection.  Even  though  several  of  their  assumptions  are  unfulfilled  in  our  case,  the 
resulting  effects  on  ion  collection  are  probably  small  enough  for  our  purposes.  We 
therefore  multiply  Eq.  (3.1)  by  the  same  factor  to  obtain  an  estimate  of  total  icn 
collection  as  influenced  by  surface-potential  effects.  The  resulting  ion-current 
dependence  on  ion  speed  ratio  is  plotted  in  Fig.  8.  For  O*  ions  at  T  =  0.1  eV 
(1160K),  H*  at  0.1  eV,  O’at  IkeV,  and  H*  at  IkeV,  we  have  S.  =  7.31,  1.83, 

0.0731,  and  0.0183  (the  latter  two  are  effectively  zero),  respectively.  The  corres¬ 
ponding  i on-current  enhancement  factors  (values  of  i.)  from  Fig.  8  are  9.50,  2.69, 
1.45,  and  1.45,  respectively. 

If  the  ambient  ions  are  H‘,  the  ion  collected  current  is  now  given  by: 

_  r  ^T.  ,  * 


Il  =  4»rs*enic{-Zi^-  }  (2.69) 

+  4,rrsJenih{  TiST  } 


where  the  subscripts  ic  and  ih  refer  to  the  cold  and  hot  ion  populations.  If  the  ions 
are  O*.  then  the  factor  2.69  in  (3.2)  should  be  replaced  by  9.50. 

The  electron  collected  current  is: 

f  i  f  1 


I.  =  4rrr 


K  ec  /^s  \ 

trsr)  “Kitt- ) 


♦  “""V  -Sr  }  '»p{  TT~  }  ' 


If  high-voltage  charging  occurs,  then  -e$  »  kT  and  the  first  term  on  the  risht- 
nana  si ae  or  tms  equation  oecomes  negligible. 


For  current  balance,  I.  =  ,re.  This  leads  to: 


2.69n 


1 .45nf^  \r 


vm./m 
i  e 


e  ■'•♦s,/kTeh 


m./me  =  43  for  H*  ions.  Therefore: 


v-y  s. 


for  H*  Ions,  with  43  and  2.69  replaced  by  172  and  9.50  for  O'  Ions.  This  is 
equivalent  to: 


el^lAT^  =  In 


hot-electron  ambient  flux 


2.69  (cold-ion  ambient  flux)  +  1.45  (hot-ion  ambient  flux) 


For  O’/H*  mixtures  and  for  hot-ion  temperatures  other  than  1  keV,  generalization 
of  this  result  is  straightforward.  Since  any  hot  ions  are  likely  to  have  W 104’ 

the  hot-ion  ambient  flux  will  exceed  the  cold-ion  ambient  flux  if  the  hot  ions 
constitute  more  than  about  IX  of  the  total  ambient-ion  number  density.  Equation 
(3.7)  indicates  that  the  onset  of  high-voltage  charging  can  be  expected  to  depend 
primarily  on  the  ratio  of  hot-electron  ambient  flux  to  the  ambient  flux  of  aU  ions,  as 
mentioned  at  the  beginning  of  this  Section.  This  completes  our  argument  in  support 
of  this  conclusion. 


In  analyzing  spacecraft  data,  one  is  therefore  likely  to  find  better  correlation  of 
spacecraft  voltages  with  the  ratio  which  appears  in  Eq.  (2.3),  or  something  nearly 
equal  to  it,  than  with  any  other  measurable  quantity,  such  as  electron  or  ion  density 
or  average  energy,  taken  individually.  This  expectation  has  been  borne  out  in  recent 
work  by  Gussennoven  et  al  (1985,  Fig.  7),  involving  charging  data  from  the  DMSP 
F6  and  F7  satelhtes.  In  calculating  values  of  this  ratio,  the  ambient  fluxes  which 
are  involved  need  to  have  been  measured  simultaneously  on  the  same  spacecraft. 

Even  though  the  approximations  made  in  deriving  (3.7)  are  severe,  an a  the  precise 
dependence  of  spacecraft  voltages  on  this  ratio  may  therefore  differ  substantially 
from  that  given  in  Eq.  (3.7)  (arid  the  coefficients  in  (3.7)  will  need  to  be  modified  if 
O*  dominates),  our  general  conclusion,  i.e.  that  spacecraft  voltages  should  correlate 
most  strongly  with  this  ratio,  or  something  nearly  equal  to  it,  is  likely  to  remain 
valid.  Furthermore,  the  dependence  of  spacecraft  voltages  on  this  flux  ratio  is  likely 
to  retain  an  approximately  exponential  form.  In  situations  where  most  secondary  and 
backscattered  electrons  emitted  by  the  spacecraft  will  escape  (see  Sec.  2),  primary- 
electron  incident  fluxes  will  be  approximately  cancelled  for  many  spacecraft 
materials  by  electron  escape  at  incident  energies  14)  to  a  few  keV  (Laframboise  et 
al,  1982a,b;  Laframboise  and  Kamitsuma,  1983;  Lai  et  al,  1983),  so  the  hot- 
electron  ambient  flux  term  in  (3.7)  needs  to  be  modified  accordingly. 

The  most  serious  approximation  made  in  deriving  (3.7)  is  probably  item  (4)  in  the 
list  at  the  beginning  or  this  Section.  This  is  because  ion  fluxes  on  downstream 
surfaces  are  likely  to  be  very  much  smaller  than  their  average  over  the  entire 
spacecraft.  They  are  also  likely  to  be  strongly  dependent  on  spacecraft  geometry, 
local  surface  potential  distribution,  and  OVH*  concentration  ratio.  Therefore,  the 
critical  value  of  ambient  flux  ratio,  at  which  the  onset  of  high-voitage  charging 
occurs,  is  likely  to  vary  substantially  among  spacecraft  having  different  geometries 
and  surface  materials.  In  particular,  for  spacecraft  having  electrically-isolated 
downstream  surfaces,  this  critical  ratio  is  likely,  because  of  local  charging  on  these 
surfaces,  to  be  much  lower  than  for  spacecraft  which  have  an  entirely  conductive 
surface  (Sec.  4). 

Furthermore,  in  contrast  with  tne  situation  for  total  ion  collection,  there  is  no 
known,  simole,  reliable  method  for  estimating  ion  fluxes  on  downstream  surfaces. 
Parks  and  Katz  (1983a,b)  have  develooed  an  ion  flux  calculation  for  the  downstream 
point  on  a  sphere  in  a  potential  which  has  a  given,  simple  analytic  form.  Detailed 
numerical  simulation,  wmch  includes  realistic  self-consistent  spacecraft  sheath 
potential  distributions,  and  which  probably  needs  to  involve  at  least  some  ion  orbit- 
foilowtng,  therefore  appears  to  be  essential.  In  Sec.  4,  we  report  on  preliminary 
results  from  a  calculation  of  this  kind. 
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So  far,  we  have  not  mentioned  the  difficulties  which  can  arise  in  measuring  the 
ambient  ion  fluxes  which  appear  in  Eq.  (3.7).  So  far,  we  have  also  defined  "ambient 
flux"  to  be  that  measured  in  an  Earth-fixed  reference  frame.  The  alternative  would 
be  to  define  it  as  that  measured  in  the  spacecraft  frame,  l.e.,  including  ram  effects. 
Ion  fluxes  measured  by  spacecraft  instruments  are  strongly  Influenced  t>y  ram 
effects.  In  fact,  the  numerical  factors  2.69,  1.45,  and  9.50,  which  appear  In  Eq. 

(3.7)  and  the  associated  discussion,  already  constitute  a  rough  ram-effect  correction, 
but  for  total  current  to  a  sphere,  not  for  local  collection  by  a  forward-facing  instru¬ 
ment  aperture.  !t  may  happen  that  the  ram -effect  correction  factors  for  an 
instrument  are  nearly  equal  to  the  above  factors,  so  that  the  instrument  measure¬ 
ment,  without  any  correction,  already  gives  a  good  estimate  of  the  denominator  of  Eq. 

(3.7) .  In  any  case,  the  response  of  the  instrument  will  depend  on  its  geometry,  and 
this  problem  has  already  been  treated  by  other  authors  (Parker,  1970;  Parker  and 
Whipple,  1970;  Whipple  et  al,  1974;  Chang  et  al,  1979;  Singh  and  Baugher,  1981; 
Comfort  et  al,  1982;  Laframboise,  1983b),  so  we  do  not  discuss  it  here. 


Parks  and  Katz  (1981)  and  Katz  and  Parks  (1983)  have  estimated  charging 

Sitentials  on  spherical  spacecraft  of  0.5m  and  5m  radius,  assuming  that  the  ions  are 
\  the  hot  electron  temperature  is  5  keV,  and  spacecraft  speed  is  8  km/sec. 

Their  results  can  be  compared  directly  with  those  given  by  our  Eqs.  (3.5)  -  (3.7). 
They  have  used  the  theory  of  Langmuir  and  Blodgett  (1924)  to  obtain  values  for 
sheath  radius  as  a  ftrction  of  spacecraft  potential.  They  present  spacecraft 
potentials  as  functions  of  the  ratio  k  of  hot  ("precipitating")  electron  ram  current  to 
ion  ram  current.  Tc  maxe  a  comparison,  their  value  of  *■  needs  to  be  expressed  in 
terms  of  our  ambient  flux  ratio.  They  have  assumed  the  ambient  electron  flux  to  be 
unidirectional.  To  convert  to  an  equivalent  isotropic  flux,  we  note  that  current  to  a 


sphere  =  4*r 


x  Isotropic  (random)  flux,  but  =  wr^2 


- r._  , - - - ,  —  ,  x  tndirectional  (ram)  flux. 

3  5 

Therefore,  equivalent  isotropic  flux  =  J  x  ltu directional  flux,  for  a  sphere. 


Also  for  a  sphere,  the  ratio  of  Ion  ram  to  random  currents  is 
U/(8kT./irm.)^  =  i  V«r  S..  Using  S4  =  7.31,  this  ratio  =  6.48,  so  therefore: 


their  k  - 


hot  electron  ram  current 

57W 


irr 


x  total  ion  random  current 
x  hot  electron  ram  flux 
T 


2 
S 

6T48'x"3'  nr  7  x  total  ion  random  flux 


(3.8) 


5778 


hot  electron  (equivalent)  random  flux 
total  ion  random  flux 


-  g  x  our  flux  ratio  R. 

With  coefficients  for  O*  used,  our  Eq.  (3.6)  gives: 

4>s  -  -50C0  in  (R/9.50).  (3.9) 

Figure  9  shows  our  result  and  theirs  [from  their  Fig.  3  (1981)  or  Fig.  2  (  1983)], 
plotted  together.  At  larger  potentials,  the  combined  set  of  results  shows  a  monotonic 
progression  toward  increased  charging  for  larger  spacecraft.  For  -<Ps  <  350V,  their 

5m  sphere  shows  more  charging  than  our  large-radius- limit  sphere.  This  is  because 
their  ion-current  enhancement  factor,  which  ;s  determined  by  the  size  of  a  sharp- 
edged  Langmuir'-Blodeett  sheath,  falls  below  ours,  which  includes  the  effect  of  a 
quasineutral  presheath.  This  discussion  suggests  that  the  tendency  toward  high- 
voltage  charging  always  increases  with  spacecraft  size,  but  magnetic-field  effects 
may  modify  this  (Laframboise,  1983a,  Sec.  1).  The  corresponding  curves  for  local 
charging,  on  surfaces  in  a  spacecraft  wake,  will  lie  to  the  left  of  those  shown  lnTIg. 
9,  but  these  remain  to  be  computed  numerically,  as  we  have  done  for  a  particular 
spacecraft  geometry  in  Sec.  4. 


4.  RESULTS  FROM  A  NUMERICAL  SIMULATION 


The  simple  treatment  developed  In  Sec.  3  led  to  a  suggestion  that  the  onset  of  high 
voltage  charging  in  low  polar  ort>it  can  be  expected  to  depend  primarily  on  the  ratio  R 
of  the  hot -electron  ambient  fltet  to  the  ambient  flux  of  all  ions.  In  this  Section  we 
present  preliminary  results  from  numerical  simulation  work  which  Is  directed 
toward  verifying  this  suggestion.  This  work  uses  an  adaptation  of  a  previously- 
existing  simulation  program  written  by  one  of  us  (Parker,  1983)  for  an  infinite- 
cylindrical  spacecraft  geometry  in  a  collisionless  plasma  crossflow.  If  the 
spacecraft  surface  is  entirely  dielectric,  then  at  progressively  Increasing  values  of 
R,  one  might  expect  surfaces  in  the  spacecraft  wake  to  be  the  first  to  undergo  high- 
voltage  charging,  followed  by  sideways-facing  and  finally  frontward-facing  surfaces, 
as  shown  schematically  in  Fig.  10.  The  numerical  results  show  that  at  least  one 
implied  feature  of  Fig.  10,  namely  the  apparent  monotonic  progression  to  larger 
negative  voltages  as  one  moves  arotnd  the  surface  of  the  cylinder  toward  the  wake 
point,  is  wrong;  we  will  return  to  this  question  later. 

Our  simulation  geometry  is  shown  In  Fig.  1 1.  Our  original  intention  was  to 
model  a  completely  dielectric  cylinder,  but  we  have  substituted  a  set  of  electrically- 
isolated  conductive  sectors  in  the  wake  region  because  we  found  that  strong  ion 
focusing  effects  occurred  in  the  wake,  and  these  produced  very  localized  ion-current 
deposition  regions  or  "ion  hot  spots",  whose  location  was  very  sensitive  to  small 
variations  which  occurred  in  the  sheath  potential  distribution  as  iteration  proceeded 
toward  a  self-consistent  set  of  surface  potentials.  Averaging  these  ion  currents  over 
the  10°  or  20°  intervals  shown  in  Fig.  1 1  suporessed  the  resulting  instabilities  and 
allowed  the  iterative  procedure  to  converge.  The  parameter  values  chosen  for  this 
study  included:  Debye  length  XQ/spacecraft  radius  rg  =0.001,  Ion  speed  ratio  S; 

(  =  drift  speed  U/ion  most-probable  thermal  speed  VZkT/mp  =  8,  ions  O*  at  a 

temperature  of  0.2  eV,  cold  electrons  at  0.2  eV,  and  hot  (auroral)  electrons  at  5 
keV.  We  have  again  made  the  assumptions  (1),  (2),  and  (3)  listed  in  Sec.  3. 

The  discretization  used  for  position  space  is  indicated  in  Fig.  1 1 .  "Inslde-out"  ion 
orbit-following  was  used  for  calculating  ion  surface  current  densities.  The  ion 
velocity-space  discretization  used  was  similar  to  that  described  by  Parker  (1977). 

To  achieve  sufficient  accuracy  in  ion  flux  values  at  wake-side  surface  points  required 
use  of  1024  incident  ion  directions  at  each  of  32  ion  energy  levels  at  each  such 
point.  The  iteration  was  started  using  an  ion  density  distribution  based  on  assuming 
that  ions  behaved  as  neutrals,  and  iterating  until  the  surface  potential,  surface  ion 
current  density,  sheath  potential  and  sheath  electron  density  distributions  all 
converged.  Electron  densities  and  currents  were  described  by  superpositions  of 
Boltzmann  factors.  Each  such  calculation  took  about  6  hours  on  the  AFGL  Cyber 
850.  The  next  intended  step  was  to  calculate  an  ion  density  distribution  corres¬ 
ponding  to  the  resulting  sheath  potential  distribution,  then  ’freeze"  this  ion  density 
distribution  and  iterate  the  other  quantities  involved  as  before,  and  so  on.  Each 
such  step  would  take  about  24  hours  on  the  same  computer,  so  we  have  not  so  far 
made  such  calculations,  and  the  results  presented  here  are  therefore  based  on  the 
"neutral  approximation"  for  ion  densities,  but  on  ion  orbit-following  for  ion  current 
collection  on  surfaces.  Further  details  of  the  calculation  method  are  to  be  published 
in  a  later  paper  (Laframboise  and  Parker,  to  be  published). 

Preliminary  results  from  these  calculations  are  shown  in  Fig.  12,  for  flux  ratios 
R  =  0. 1,  0.25,  0.5,  and  5.0.  On  the  front  and  sides  of  the  cylinder  (0  S  100°), 
surface  potentials  seem  to  increase  more-or-less  monotonicaily  as  one  moves  away 
from  the  front  (0  =  0°),  except  for  some  relatively  small  oscillations  which  are 
probably  spurious,  and  may  be  caused  either  by  two  coarse  a  position-space 
discretization,  incomplete  convergence,  or  insufficiently  fine  sampling  of  the  ion 
velocity  space  (even  with  32,768  orbits  followed  per  surface  point,  very  few  of  these 
will  connect  back  to  the  ambient  plasma  with  an  ambient  velocity  close  to  that  of  the 
heavily-populated  part  of  the  ion  distribution). 

However,  on  rearward  surfaces  (0^1 10°),  there  is  some  non-monotonicity  which 
appears  to  be  reai.  For  R  =0.1,  0.25,  and  0.5,  there  are  three  surface-potential 
maxima,  one  at  the  rearmost  point  (0  =  180°)  and  one  on  each  side  of  it.  For  R  = 
0.1  and  0.25,  these  features  were  almost  unchanged  (the  potentials  of  the  maxima 
changed  by  less  than  2%)  when  as  few  as  4,096  ion  orbits  (4  ion  energy  levels)  were 
usea,  and  this  attests  to  their  reality.  For  R  =  5,  we  obtain  four  maxima,  two  on 
either  sice  ot  the  rearmost  point.  The  minima  between  these  Deans  correspond  to 
the  ion  'hot  spots  ,  or  deposition  points  of  highly  beam-like  ion  populations,  men¬ 
tioned  above.  Most  real  situations  would  not  possess  the  symmetry  about  the  fore-aft 
line  which  our  orobiem  does,  but  our  results  nonetheless  indicate  that  high-voitage 
wake  regions  of  spacecraft  are  likely  to  contain  very  oeam-like  ion  population 
components  generally. 
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Trus  is  related  also  to  the  fact  that  the  ambient  ion  distribution  is  a  highly-direc¬ 
ted  (Sj  =  8)  one.  A  situation  with  ion  drift  is  very  different  than  one  without.  A 

drifting  distribution  is  not  an  equilibrium  one.  The  drift  provides  the  ions  with 
thermodynamic  free  energy  which  can  support  a  much  greater  variety  of  phenomena 
than  in  the  noodrifting  case,  including  self-excited  oscillations  (Kr-all  and 
Trivelpiece,  1973,  Crs.  2  and  9).  For  related  reasons,  the  drift  also  makes  possible 
ion  focusing  effects  including  those  discussed  here.  In  the  steady-state  treatment 
described  here,  self-excited  oscillations  cannot  be  treated.  However,  oscillations 
recently  observed  in  the  outer  portions  of  the  Shuttle  wake  involve  density  fluctu¬ 
ations  of  only  a  few  per  cent  (Murphy,  1985),  so  they  probably  do  not  have  much 
effect  on  the  steady~state  wake  properties  studied  here. 

Beam-like  ion  components  have  important  implications  for  more  complicated 
situations,  such  as  the  interior  of  the  Shuttle’s  cargo  bay.  As  our  work  has 
indicated,  beam  impact  points  will  be  very  sensitive  to  details  of  the  potential  distri¬ 
bution  on  or  near  a  spacecraft.  Ions  may  come  around  a  comer  of  the  spacecraft, 
just  miss  one  equipment  box,  and  impact  on  one  small  corner  of  the  next  box,  which 
may  then  come  to  a  potential  very  different  than  these  of  its  surrourdings. 

In  addition,  more  than  one  ion  beam  may  be  present  at  some  points  in  the  sheath. 
Beams  may  come  around  comers  from  opposite  sides  of  a  spacecraft,  and  beams  may 
also  be  present  which  have  circled  the  spacecraft  one  or  more  times.  Ion 
distributions  as  a  function  of  angle  may  therefore  contain  many  narrow  peaks"  and 
"valleys".  Simulation  of  a  very  detailed  kind,  demanding  the  maximum  available 
computer  power,  will  be  necessary  to  resolve  such  situations. 

In  Fig.  9,  we  have  also  plotted  the  largest  value  of  negative  surface  potential  for 
each  of  tne  four  values  of  &  shown  in  Fig.  12,  for  compartsion  with  our  simple 
prediction  from  Sec.  3.  These  values  are  shown  as  circled  points  in  this  Figure. 
These  results  Qualitatively  resemble  the  theoretical  curves  in  the  same  Figure,  and 
also  our  expected  behaviour  shown  in  Fig.  10,  in  that  they  appear  to  show  almost  no 
dependence  on  R  until  a  threshold"  value  is  reached,  and  then  tney  show  a  rapid 
increase.  Furthermore,  this  increase  aooear-s  to  begin  at  a  substantially  lower  R 
value  than  for  these  curves,  as  one  expects  since  the  curves  are  for  whole-body 
charging  and  the  plotted  points  are  for  charging  of  electrically-isolated  surfaces  in 
the  wake. 

One  feature  of  these  results  is  unexpected:  the  "pre-threshold"  wake  charging, 
instead  of  being  close  to  zero,  is  already  several  hundred  volts.  However,  this 
apparent  thresnoid  is  not  tne  true  one.  'High-voltage  charging  conventionally  means 
that  surface  potential  exceeds  100V  in  magnitude.  For  our  conditions,  ana  for  a 
surface  potential  of  -lOCV,  the  hot -electron  flux  exceeds  the  cold  electron  flux  when 
R  >  1.2  x  1C'J;S!  Also,  for  R  =  0.1  (the  left  margin  of  Fig.  9),  if  we  calculate 
downstream-oemt  potential  using  the  neutral  approximation  for  ion  flux  [Tsien,  1946; 
Parks  and  Katz,  1983b,  Eq.  4j,  we  obtain  <*>s  =  -3.3  x  10*V.  For  R  =  0  (cold 

electrons  only),  a  similar  calculation  yields  <pg  -  -14.8V.  It  is  evident  that  all  the 

numerical  simulation  results  shown  in  Fig.  9  are  post-threshold"  in  the  sense  that 
they  involve  a  current  balance  primarily  between  the  hot  electron  fiux  and  wake- 
region  Ion  fluxes  that  are  already  heavily  modified  by  orbit  curvature  in  strong  wake- 
region  electric  fields.  This  is  evidently  a  situation  m  which  even  a  small  amount  of 
ambient  H*  can  be  expected  to  produce  a  large  decrease  in  wake-region  potentials,  and 
we  intend  to  include  H*  effects  in  future  calculations. 
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TABLE  1.  Low-Earth-Orbit  Conditions** 


(a) 

CHARACTERISTIC  LENGTHS 

Ambient  Debye  Length: 

£  1  cm 

Thickness  of  lkV  and  5kV  spherical ly-symmetric 
sheath***  arxxnd  a  sphere  of  radius  3m: 

2.5m  and  5.8m 

Thickness  of  lkV  and  5k V  planar  Child-Langmuir 
sheath: 

5.0m  and  17m 
(WHY  SO  LARGE?) 

Particle  gyroradii: 

ambient  electrons  (0.1  eV) 

2  cm 

secondary  electrons  (3  eV) 

13  cm 

auroral  electrons  (10  keV) 

8  m 

ions  (0*;  0.1  eV) 

3m  in  "rest"  frame 

27  m*  in  spacecraft  frame 
(larger  inside  a  high-voltage  sheath) 

(b) 

CHARACTERISTIC  SPEEDS: 

Ion  thermal  speed: 

Spacecraft  speed: 

1  km/sec 

£  8  km/sec, 
depending  on  orbit 

SITUATION 

>  B 

Electron  thermal  speed: 

100  km/sec, 

for  ambient  electrons 

-J 

•  KCSOTHERMAL" 

*  Particle  motions  do  not  depend  on  the  frame  of  reference  in  which  they  are  viewed. 
The  transformation  from  rest  frame  to  Spacecraft  frame  produces  a  V  x  B  electric 
field  S  0.35V/m  where  V  =  spacecraft  velocity  and  B  =  magnetic  induction;  electric 
fields  inside  a  lkV  sheath  are~  200  V/m. 

**  deLeeuw  (1967),  p.  1564;  Martin  (1974). 


***  Al’pert  et  al  (1965),  Fig.  72. 


Figure  1.  Effect  of  surface  orientation  on 
escape  of  emitted  electrons.  In  (a),  the 
spacecraft  surface  is  perpendicular  to 
tne  magnetic  field  B,  and  the  emitted 
electrons,  which  experience  an  electric 
force  -eE  directed  away  from  the  surface, 
all  escape.  In  (b),  the  spacecraft  surface 
Is  nearly  parallel  to  B,  and  almost  all 
of  the  emitted  electrons  return  to  the 
surface,  even  though  they  still  experience 
an  electric  force  directed  away  from  it. 
Note  that  the  component  of  E  perpendi¬ 
cular  to  B  resuits  only  in  an  £  x  B  drift 
parallel  to  the  surface. 
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Flpr*  8.  Dependence  of  ion  cu-rent 
to  a  sphere  on  ion  speed  ratio. 


Figure  9.  Dependence  of  spacecraft  surface  potential  on  hot  electron/total  ion 
ambient  flux  ratio.  The  four  circled  points  are  the  largest  values  of  surface  potential 
from  the  jmerical  solutions  discussed  in  Sec.  4  and  presented  in  Fig.  12. 


Figtre  10.  Expected  genera.  appearance  of  graphs  of  spacecraft  surface  potential  as 
functions  of  ambient  hot  electron  to  total  ion  flux  ratio  R. 
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EXTRA  GRID  LINES 
AT  M*  AND  OS* 


Figure  il.  Geometry  of  simulated  spacecraft  used  for  computational  purposes.  The 
outer  boundary  of  the  computational  domain  is  a  square  whose  sides  are  at 
*  =  *5rs,  y  =  *5rs,  where  rg  is  the  spacecraft  radius.  The  domain  is  covered  by  a 

finite-element  grid  in  which  there  are  10  intervals  along  each  radial  line  from  the 
spacecraft  surface  to  the  outer  boundary,  and  the  size  of  these  intervals  is  propor¬ 
tional  to  the  radius.  Within  the  gaps  between  the  conductive  sectors,  the  potential  is 
assumed  to  vary  linearly  with  surface  position  on  the  cylinder.  The  dots  on  the 
conductive  sectors  indicate  locations  on  them  where  ion  currents  are  calculated. 


Figure  12.  Surface  potentials  as  functions 
of  position  for  the  spacecraft  geometry  shown 
in  Fig.  11,  for  ratios  R  of  hot  electron  to 
total  ion  ambient  flux  as  shown. 
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CALCULATION  OF  SECXXvDARY- ELECTRON  ESCAPE  CURRENTS 
FROM  FEGATIVELY-CHARGED  SPACECRAFT  SURFACES  IN  A  MAGNETIC 

FIELD 


J.  G.  Laframboise 

Physics  Department,  York  University 
Toronto.  Canada  M3J  1P3 


ABSTRACT 


In  low  Earth  orbit,  the  geomagnetic  field  B  is  strong  enough  that  secondary 
electrons  emitted  from  spacecraft  surfaces  have  an  average  gyrcradius  much 
smaller  than  typical  dimensions  of  large  spacecraft.  This  implies  that  escape 
of  secondaries  will  be  strongly  inhibited  on  surfaces  which  are  nearly  parallel 
to  B  ,  even  if  a  repelling  electric  field  exists  outside  them.  This  effect  is 
likely  to  make  an  important  contribution  to  the  current  balance  and  hence  the 
equilibrium  potential  of  such  surfaces,  making  high-voltage  charging  of  them 
mere  likely.  We  present  numerically-calculated  escaping  secondary  electron 
fluxes  for  these  conditions.  For  use  in  numerical  spacecraft -charging 
simulations,  we  also  present  an  analytic  curve-fit  to  the  results  for  the 
important  case  of  normal  electric  field  (uniformly-charged  surfaces).  This 
curve-fit  is  accurate  to  within  3%  of  the  emitted  current.  For  strong  normal 
electric  fields,  escape  is  effectively  suppressed  only  when  a  surface  is  parallel 
to  B  within  a  few  degrees  or  less,  and  this  leads  to  "sensitivity  effects"  in 
attempts  to  predict  auroral-zone  spacecraft  charging.  A  nonzero  tangential 
component  in  the  surface  electric  field  can  greatly  enlarge  the  range  of  surface 
orientations  for  which  escape  is  suppressed,  and  can  also  produce  large  surface 
currents. 
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i.  INTRODUCTION 


The  prediction  of  high-voltage  charging  or  other  environmental  effects  on  a 
spacecraft  in  low  Earth  orbit  appears  likely  to  be  more  complicated  than  in 

geostationary  orbit,  for  at  least  three  reasons.  j 

( 

These  reasons  are:  (a)  space  charge  effects  (on  sheath  and  wake  potentials) 
are  more  important,  because  space-charge  densities  are  much  higher  (the 
Debye  length  is  no  longer  >>  typical  spacecraft  dimensions)  (b)  ion  flow 
effects  are  more  important,  because  spacecraft  orbital  speed  >  ion  thermal 
speeds  (c)  the  geomagnetic  field  B  is  likely  to  have  an  important  influence  on 
charged-particle  motions  because  B  is  now  much  larger,  and  not  all  of  the 
average  particle  gyrorariii  of  importance  are  any  longer  >>  typical  spacecraft 
dimensions. 

We  wish  to  investigate  an  important  consequence  of  (c),  which  concerns  the 
escape  of  secondary  electrons  emitted  from  spacecraft  surfaces.  Our  discus¬ 
sion  will  aiso  apply,  with  minor  modifications,  to  photoelectron  or 
backscattered-electron  escape.  In  low  Earth  orbit,  in  the  auroral-zone 
geomagnetic  field  (|B|=  0.44  gauss  =  4.4  x  1 0_5T) ,  the  gyroradius  of  a 
"typical"  3eV  secondary  electron  and  a  10  keV  auroral  electron  are  13  cm  and 
S  n,  respectively.  The  average  gyroradius  of  "cold"  ionospheric  electrons 
(temperature  T  =  0. 1  eV)  in  the  same  B  is  even  smaller  (2  cm),  but  this  is 
not  an  important  parameter  in  most  cases  because  these  electrons  are  repelled 
:f  the  spacecraft  potential  is  negative,  and  their  density  is  then  well- 
apprcximatec  by  a  3oitzmanr.  factor,  which  is  unaltered  by  B  effects. 


The  reasor.  why  B  affects  secondary-electron  escape  is  shown  in  Fig.  1.  In 
Fig.  1(a),  the  spacecraft  surface  is  perpendicular  to  B,  and  the  emitted 
electrons,  which  experience  an  electric  force  -eE  directed  away  from  the 
surface,  all  escape,  helping  to  discharge  it.  In  Fig.  1(b),  the  SDacecraft 
surface  is  nearly  parallel  to  B,  and  almost  all  of  the  emitted  electrons  return 
to  it,  even  though  they  still  experience  an  electric  force  directed  away  from  it 
Tr.e se  electrons  therefore  are  unable  to  helD  discharge  it,  so  a  surface  nearly 
parallel  to  B  is  mere  iikeiy  to  charge  to  a  large  negative  voltage.  Note  that 
the  component  of  E  which  is  perpendicular  to  B  results  only  in  an  E  x  B  drift 
parallel  to  the  surface. 

For  any  object  much  larger  than  13  cm,  the  escape  of  secondary  electrons 
will  be  strongly  affected  by  this  process.  For  example,  most  surfaces  on  the 
Shuttle  are  effectively  "infinite  planes '  by  this  criterion.  On  the  other  hand, 
the  average  gyroradius  of  high-energy  auroral  electrons  is  comoarable  to 
Shuttle  dimensions,  so  the  deposition  of  these  electrons  onto  Shuttle  surfaces 
is  likely  to  be  only  moderately  inhibited. 

For  a  larger  object  (size  >>  8  m),  deposition  of  auroral  electrons  will 
also  become  strongly  orientation-dependent,  with  both  collection  and  escape  of 
electrons  now  being  inhibited  on  surfaces  nearly  parallel  to  B.  This  suggests 
that  high-voltage  charging  of  such  surfaces  may  be  mere  likely  on  ob-ects  of 
intermediate  size  than  on  either  larger  or  smaller  ones.  In  tke  calculation  of 

Parks  and  Katz  (1981),  Katz  and  Parks  (1983),  the  tendency  toward 
nigh-voltage  charging  increased  with  spacecraft  size  because  in  their  model. 


ion  collection  increased  less  rapidly  with  spacecraft  size  than  did  electron 
collection.  To  determine  which  of  these  two  effects  predominates  will  require 
more  detailed  calculations  than  have  been  done  so  far. 

As  already  mentioned,  strong  ion  flow  effects  also  are  generally  present  in 
low  orbit;  the  ion  speed  ratios  (flow  speed/most  probable  ion  thermal  speed) 
for  hi'*'  at  1  keV,  at  0.1  eV,  and  0+  at  0.1  eV  are  0.02,  1.8,  and  7.3, 
resoectiveiy.  Whenever  the  latter  is  the  predominant  ion  species,  ion 
collection  on  downstream  surfaces  will  therefore  be  strongly  inhibited.  If  a 
surface  is  simultaneously  downstream  and  nearly  parallel  to  B,  as  is  likely  to 
be  the  case  in  the  auroral  zones,  then  the  tendency  for  high-voltage  charging  to 
occur  on  it  will  be  greatly  increased  (Fig.  2). 

To  "straightforwardly"  include  B  effects  on  secondary  electron  emission  in  a 
large  two  or  three  dimensional  simulation  program  would  involve  the 
numerical  integration  of  very  large  numbers  of  secondary-electron  orbits.  The 
resulting  computing  costs  usually  would  be  formidable,  especially  since  these 
orbits  would  have  relatively  large  curvatures.  A  desirable  alternative  is  to 
"parameterize"  the  situation  by  treating  in  advance  a  simplified  but  still 
sufficiently  realistic  model  problem.  In  order  to  do  this,  we  make  the  approx¬ 
imations  described  in  the  next  Section. 


2.  THEORY  FOR  E  NORMAL  TO  SURFACE 


We  assume  that  the  spacecraft  surface  is  an  infinite  plane,  and  the  electric 
and  magnetic  fields  E  and  B  outside  it  are  uniform.  In.  this  Section,  we  aiso 
assume  that  the  electric  force  -eE  on  electrons  is  directed  along  the  outward 
normal  to  the  surface;  here  e  is  the  magnitude  of  the  elementary  charge.  This 
assumption  is  relaxed  in  Sec.  4,  in  order  to  permit  variations  of  potential 
along  the  surface  to  be  taken  into  account.  We  assume  that  the  secondary 
electrons  are  emitted  with  a  Maxwellian  distribution  corresponding  to  a 
temperature  T.  The  ratio  i  =  I/I  of  escaping  to  emitted  flux  is  then  a  funct¬ 
ion  of  two  parameters:  the  angle  0  between  the  surface  normal  and  the  direct¬ 
ion  of  B  (Fig.  3),  and  a  parameter  describing  the  strength  of  E.  A  convenient 
choice  for  this  parameter  is  the  difference  in  potential  across  a  mean  second¬ 
ary-electron  gyroradius  a  =  (i/eB)(7rmkT/2)^,  divided  by  kT/e,  where  m  is 
electron  mass  and  k  is  Boltzmann’s  constant. 


This  quotient  is: 


where  E  =  |  E  |  and  B  =  |  B  | . 


This  quantity  aiso  has  an  alternative,  more  useful  interpretation:  it  is  the 
ratio  of  the  magnitude  |E  x  B  \/B2  of  the  E  x  B  drift  speed,  to  one-half  the 


2  3 


SfV 


mean  thermal  speed  (8kT/7rm)”  of  the  emitted  electrons.  It  is  useful  to 
estimate  the  value  of  e  for  a  high-voltage  spacecraft  sheath  in  low-orbit 
conditions.  To  do  this,  we  use  the  sheath  solution  of  Al’pert  et  al  (1965, 

Table  XXIV  and  Fig.  72).  For  a  1  kV  and  a  5  kV  sheath  around  a  sphere  of 
radius  3m  in  a  ccllisicr.iess  plasma  having  an  ambient  ion  temperature  of 
C.l  eV,  number  density  of  3  x  10s  cm'3,  and  resultant  (ion)  Debye  length  of 
3.43  cm,  their  resuits  give,  respectively,  sheath  thicknesses  of  2.6  and  6.1  m, 
and  surface  electric  fields  E  =  0.86  and  2.9  kV/m.  Using  B  =  4.4  x  10'5T 
and  T  =  3  eV  for  secondary  electrons,  we  then  obtain  e  =  33.9  and  1 14.2. 

Beth  of  these  are  relatively  large  values,  whose  significance  can  be  understood 
if  we  consider  what  would  happen  if  e  were  infinite. 

In  this  limit,  it  is  easy  to  show  that  secondary  electrons  would  all  escape 

unless  B  were  exactly  parallel  to  the  surface  (6  were  90°).  This  can  be 
shown  as  follows.  In  this  limit,  secondary  electrons  would  have  no  "thermal" 
motion.  The  (y,z)  projection  of  their  motion  would  then  be  similar  to  that 
shewn  in  Fig.  4.  This  motion  would  be  the  sum  of:  (i)  an  E  x  B  drift  in  the 
y  direction  (ii)  a  uniform  acceleration  along  B  ,  whose  projection  in  the  (y,z) 
plane  would  be  upward  (iii)  just  enough  gyromotion  to  produce  a  cycloidal  path 
when  combined  with  (i),  so  that  in  the  absence  of  (ii),  the  electron  would  (just) 
return  to  the  surface  at  the  end  of  each  gyroperiod.  In  the  presence  of  (ii), 
these  "return  points"  are  displaced  upward  by  progressively  increasing  amounts 
(Fig.  4),  so  the  electron  can  never  return  to  the  surface,  unless  B  is  exactly 
parallel  to  the  surface,  so  that  the  upward  component  of  -eE  along  B  vanishes. 

This  result  suggests  that  for  large  finite  values  of  e  (including  the  values 


calculated  above),  electron  escape  is  likely  to  be  almost  complete  except  for  6 
very  near  90°,  where  it  should  drop  to  zero  very  steeply.  The  occurrence  of 
high-voltage  charging  in  marginal  circumstances  may  therefore  depend  very 
strongly  on  the  precise  orientation  of  a  surface. 

The  escaping  secondary-electron  flux  is  given  by: 

I  =  III  f(v0)  H(v)  voz  d»V0 
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(2.2) 

where:  v^  is  the  initial  velocity  of  an  emitted  electron,  f(vQ)  =  d-Vi/d^  is 
the  velocity  distribution  of  emitted  electrons,  n  is  a  reference  number  density, 
and  H(vq)  is  equal  to  1  for  escaping  electrons  and  0  for  those  which  return  to 
the  surface.  The  emitted  flux  is: 

I  =  n(kT/277m)^.  (2.3) 

o 

We  also  introduce  the  dimensionless  velocity: 

u  =  v  (m/2kT)^  .  (2.4) 


Equation  (2.2)  then  becomes: 
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which  is  in  a  form  suitable  for  numerical  summation.  The  quantities  Ujim 

u,.  ,  u  i.  ■  are  the  values  of  u  for  which  H  changes  between  0 

and  1  for  each  u  and  u  .  These  values  must  be  found  by  numerically 
ox  oy  J  J 

determining  which  particle  orbits  reimpact  the  surface.  These  orbits  can, 
however,  be  determined  in  analytic  form,  with  time  as  a  parameter.  To  do 
this,  we  use  the  coordinate  system  shown  in  Fig.  3,  together  with  a  y-axis  (not 
shewn)  directed  into  the  plane  of  the  Figure.  The  equation  of  motion  for  an 


electron  is: 


v  =  -  -  (E  +  v  x  B  ). 
m 


(2.6) 


We  solve  this  with  the  initial  conditions  t  =  y  =  r>  =  0,  vt  =  v  e,  v  =  v  ,  and 

1  '  1  $  oq  y  oy 

v  =  v  .  We  introduce  the  dimensionless  variables: 
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X  -  x/a,  y  =  y/a,  etc; 


(2.7) 


r  =  u>  t  =  feB/m)t. 
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In  the  present  discussion,  e  and  e  are  both  zero,  but  for  later  use  (Secs. 

x  y 

4  and  5),  we  have  retained  these  quantities  in  the  formulas  below.  We  obtain: 


u  »  =  u  w  sin  0  +  u  cos  0; 
o£  ox  oz 
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z  =  £  cos  0  +  q  sin  0. 

•v 

Equations  (2.8)  can  also  be  differentiated  to  find  dz/dr.  The  numerical 
procedure  for  finding  the  quantities  Ujim  ^  in  Eq.  (2.5)  then  involves 

<v  % 

calculating  z  and  dz/dr  at  a  succession  of  points  along  an  orbit  (the  electron 

will  reimpact  during  the  first  gyroperiod  0  <  r  £  2rr  if  at  all,  so  this  interval 

always  suffices),  and  making  the  appropriate  tests  on  these  quantities  to  find 

out  whether  the  orbit  reimpacts  or  escapes.  For  each  uo^  .  and  u^  ,  this  is 

done  for  a  succession  of  values  of  u  .  These  tests  also  vield  the  local 

OZ  ' 

minimum  of  z(r)  if  one  exists.  Whenever  a  change  occurs  between  no  escape 
and  escape  from  one  such  value  of  uq2  to  the  next,  an  interpolation  using  these 

minima  can  be  used  to  provide  the  corresponding  value  of  u^m  In  cases 
where  they  are  unavailable,  the  arithmetic  mean  of  the  two  successive  u 
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values  is  used. 

We  have  chosen  the  abscissas  u  .  and  u  in  Eq.  (2.5)  by  first  solving 
numerically  the  transcendental  equation  $  +  $  erf(u^/V2)  =  k/n  for 

k  =  1,2, _ ,  n-1.  The  resulting  values  u^  then  subdivide  the  normalized 

Gaussian  distribution  exp(-u2)/V7r  into  n  "slices"  whose  areas  decrease  as 
exp(-4u^2)  when  n  is  large.  In  terms  of  these  u^,  we  then  choose  values 
uQ  k  at  the  centroids  of  these  slices,  and  weights  ck  equal  to  the  areas  under 
each.  We  obtain: 


exp(-uk2)  -  exp(u^ _^2) 

u  -  - 

’  Virterf  u^  -  erf 

ck  =Merfuk-erfuk.1) 


(2.9) 


(2.10) 


for  k  =  1,2, _ ,n.  A  convenient  method  for  calculating  the  required  values  of 

erf  u  =  1-erfc  u  has  been  given  by  Shepherd  and  Laframboise  (1981).  The 
resulting  vaiues  u  ,  are  then  used  to  provide  the  required  values  of  u 

0,K  OX^i 

and  u  . ,  and  the  c,  are  used  to  provide  values  of  (l/Vrr)  Au  .  exp(-u  . 2) 
oy^  j  k  “  ox  ^  1  ox  y  i 

and  (l/Vrr)Au  .  exp(-u  .2),  for  use  in  Eq.  (2.5).  We  have  provided  values 
oy,j  oy,j 

of  u^z  for  use  in  determining  the  Uj.^  k  vaiues  in  Eq.  (2.5)  by  solving  the 
equation  l-exp(-u  ,  2/2)  =  k/n  for  k  =  0,1,2,...,  n-1.  This  gives  u^  ,  = 

*  OZ,K  0Z,K 

* 

(  2in[  1/(1  -k/n) ]}J;  these  values  are  distributed  most  densely  near  uq^  =  0,  but 
still  densely  enc-g.n  at  large  uqz  that  the  resulting  intervals  give  vanishing  flux 
contributions  in  this  limit.  This  completes  the  definition  cf  the  procedure  used 
for  calculating  the  ratio  I/I  of  escaping  to  emitted  flux. 
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3.  RESULTS  AND  DISCUSSION  FOR  E  NORMAL  TO  SURFACE 


Escaping  secondary-electron  current  densities,  computed  as  described  in 
Sec.  2,  are  shown  in  Table  1  and  Fig.  5.  Each  value  of  i  -  I/I  was  calculated 
using  80x80x40  orbits,  whose  initial  velocity  components  u0xiu0y>  and  uqz 
were  chosen  as  described  in  Sec.  2,  and  with  points  on  the  orbits  calculated 
at  intervals  Ar  =  tt/45.  For  8  values  of  e  and  1  1  values  of  0,  the  resulting 
calculation  took  about  100  hr  total  on  a  Hewlett-Packard  1000F  minicomputer 
with  Vector  Instruction  Set.  The  results  are  accurate  to  within  about  0.2%  or 
better.  The  result  for  e  =  0  is  just  the  analytic  result  i  =  cos  0.  To  see  why 
this  is  so,  we  consider  the  electron  orbit  shown  in  Fig.  6,  which  has  been 
fictitiously  extended  so  as  to  pass  through  the  surface  and  re-emerge  from  it. 
In  the  absence  of  an  electric  field  (e  =  0),  this  orbit  has  the  same  speed  at  the 
re-emergence  point  C  as  at  the  emission  point  A.  Since  we  have  aiso  assumed 
that  the  emitted  velocity  distribution  is  isotropic,  and  therefore  a  function  of 
speed  only,  the  real  orbit,  for  which  C  is  the  emission  point,  must  carry  the 
same  population  as  would  the  fictitious  re-emerged  orbit.  The  flux  crossing 
the  reference  surface  DE,  which  is  1  B,  is  therefore  the  same  as  if  such 
passages  and  re-emergences  actually  occurred,  and  is  the  same  as  if  another 
reference  surface  FG,  also  i  B,  were  emitting  electrons  having  the  same 
velocity  distribution.  However,  in  reality,  the  electrons  come  from  the  real 
surface  HJ,  which  is  not  i  B,  and  all  the  electron-orbit  guiding  centers  which 
are  inside  any  given  magnetic-flux  tube  through  DE  will  also  be  inside  the 


29 


«'i 


1 


ore  section  of  the  same  flux  tube  onto  HJ,  anc  the  ratio  of  the  intersection  areas 
of  this  tube  with  HJ  and  DE  is  just  sec  0.  The  ratio  of  escaping  to  emitted 
flux  must  therefore  be  the  reciprocal  of  this,  or  cos  0,  as  stated  above. 

Also  evident  in  F:g.  5  is  the  fact,  mentioned  in  Sec.  2,  that  when  e  is  large 
enough,  electron  escape  becomes  essentially  complete  except  when  0  is  very 
nearly  90°.  This  means  that  in  Shuttle  high-voltage  charging  conditions,  for 
which  30  ^  e  £  120  (Sec.  2),  the  occurrence  of  high-voltage  charging  in 
marginal  circumstances  may  depend  very  strongly  on  the  precise  orientation  of 
a  surface.  A  slowly-rotating  surface  which  passes  through  tangentiality  to  B 
may  experience  a  sudden,  brief  high-voltage  charging  event.  For  the  same 
-eason,  attempts  to  predict  high-voltage  charging  may  be  afflicted  by 
sensitivity"  problems:  if  one  attempts  to  predict  "worst-case"  charging  by 
assuming  that  secondary  electrons  do  not  escape,  then  the  resulting  predictions 
are  iikely  to  be  overly  pessimistic  most  of  the  time.  On  the  other  hand,  if  one 
assumes  that  secondaries  do  escape,  correct  predictions  will  be  obtained 
almost  all  of  the  time,  but  occasionally  a  large  underestimate  of  charging  will 
occur. 

In  a  real  situation,  E  would  not  be  uniform,  but  would  decrease  with  distance 
from  the  surface,  contrary  to  our  assumptions.  Our  results  can  therefore  be 
expected  to  overestimate  electron  escape.  This  would  probably  not  be  a  large 
effect,  but  this  presumption  remains  to  be  verified.  An  approximate  compen¬ 
sation  for  it  can  be  made  by  calculating  e  using  an  electric  field  value  which  is 
averageci  over  the  first  mean  gyroradius  distance  from  the  surface. 
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The  results  in  Table  1  are  approximated  to  within  2.5%  of  I  by  the 
empirical  formula: 

r 

a  =  1  +  1.35e‘-‘39«  exp'  0.083725  jl+tanh  [1.9732  in  (y^)]  • 


-0.07825  in  [1  +  (e/8.5) ,-78148]  [; 


b  =  0.38033c0-95892  exp ^2.0988^1  +  tanh  [1.49  in  (3735))}  ); 

c  =  in  (9O°/0); 
i  =  cos  [90°  exp(-ac-bc2)] . 


(3.1) 


This  formula  also  has  the  correct  limiting  behavior  when  e  ->0  or  co,  or 
0  0°  or  90°.  An  approximation  formula  for  the  emitted  flux  is  also  avail¬ 

able  [Eqs.  (5)  and  (6)  of  Laframboise  et  al  (1982a),  and  Laframboise  and 
Kamitsuma  (1983)].  The  results  presented  in  this  Section  have  already  been 
presented  by  Laframboise  (1985). 

4.  THEORY  FOR  NON-NORMAL  DIRECTIONS  OF  E 

If  a  spacecraft  surface  is  charged  to  a  nonuniform  potential,  the  electric  field 
E  outside  it  has  a  nonzero  component  tangential  to  the  surface.  The  resulting 
problem  geometry  is  more  complicated,  and  two  additional  angles,  cr  and  \</, 
which  define  the  direction  of  the  electric  force  -eE,  need  to  be  specified 
(Fig.  7). 

It  may  then  happen  that  even  though  -eE  is  still  directed  away  from  the 
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surface,  its  projection  along  the  direction  of  B  is  directed  toward  the  surface, 
as  shown  in  Fig.  8a.  This  happens  when  the  angle  between  -eE  and  either  B  or 
-B,  whichever  is  outward,  is  greater  than  90°.  The  resulting  situation  is 
analogous  to  that  for  a  sailboat  tacking  against  the  wind  (Fig.  8b):  even  though 
the  resultant  wind  force  on  the  sail  has  a  downwind  component,  its  projection 
parallel  to  the  boat’s  direction  of  motion  (roughly  along  its  centreline)  has  an 
upwind  component,  and  the  boat  therefore  moves  upwind. 

When  this  situation  exists,  every  emitted  electron  returns  to  the  surface, 
and  the  escaping  flux  is  zero.  The  resulting  situation  then  differs  from  that 
shown  in  Fig.  5  in  at  least  one  important  respect.  The  range  of  surface 
orientations  (relative  to  B)  for  which  electron  escape  is  entirely  prevented,  is 
no  longer  infinitesimal  (at  Q  =  90°),  but  finite,  and  this  may  greatly  enlarge 
the  portion  of  a  spacecraft’s  surface  for  which  secondary-electron  escape  is  not 
available  as  a  discharge  mechanism. 

Even  though  no  electrons  escape  under  these  conditions,  the  possibility 
exists  that  they  may  travel  a  relatively  long  distance,  equal  to  many  gyro- 
radii,  parallel  to  the  surface  before  returning.  This  may  produce  relatively 
large  surface  currents,  and  these  may  modify  substantially  the  charge  distrib¬ 
ution  on  the  spacecraft.  The  question  of  surface  currents  resulting  from  migr¬ 
ation  of  reimpacting  electrons  along  the  surface  is  to  be  examined  in  a  sub- 
secuent  caoer  (J.G.  Laframboise.  to  be  oublisned).  Here,  we  ooint  out  onlv  one 

.  .  ;  '  ‘  '  '  i  J 

general  feature  of  this  charge  migration.  The  general  motion  of  an  emitted 


electron  is  a  superposition  of  gyromotion  about  B,  an  E  x  B  drift,  and  an 

acceleration  in  the  -eE  direction.  The  ExB  drift  direction  may  have  either  an 

inward  or  outward  normal  component.  If  it  is  outward,  migration  distances 

over  the  surface  will  be  larger,  and  so  therefore  will  surface  currents.  Since 

B  is  in  the  (x,z)  plane,  the  z-component  of  Ex  B/B2  is  -B^  E y/B2.  For 

B  >0,  this  component  has  a  sign  opposite  to  that  of  E  .  Therefore,  when 
x  y 

Ey  <  0  [region  (3)  in  Fig.  9],  the  migration  distance  of  reimpacting  elec¬ 
trons,  and  therefore  also  the  surface  currents  produced  by  them,  are  likely  to 
be  much  larger  than  when  E^  >  0  [region  (2)  in  Fig.  9].  These  surface 
currents  will  be  primarily  along  the  tangential  projection  of  the  ExB  drift 
direction,  rather  than  of  -eE,  so  they  may  have  little  effect  on  helping  to 
discharge  the  spacecraft. 

These  predictions  are  based  on  the  assumption  that  E  and  B  are  spatially 
uniform.  However,  if  average  migration  distances  become  equal  to  many 
gyroradii,  it  is  then  more  likely  that  this  assumption  will  be  seriously  in 
error,  and  many  of  our  predictions,  especially  the  quantitative  ones  presented 
in  Sec.  5,  may  then  become  unreliable. 


Before  we  present  computed  results,  we  need  to  determine,  in  terms  of  the 
angles  0,  a,  and  in  Fig.  7,  when  the  projection  of  -eE  along  B  is  directed 
toward  the  surface.  This  projection  is: 


(-eE»B)B/B2 


(4.1) 


The  z -component  of  this  is: 

-e(E  B  +  E  B  )  B  /B2 

XX  Z  Z'  z 

=  -eE  cos2©  cos  cr  (1  +  tan  a  cos  \p  tan  0)  (4.2) 

Since  -eE  cos2  0  cos  a  >  0,  this  means  that  escape  is  prevented  if: 

tan  a  cos  v//  tan  0  <  - 1 , 
of  equivalently: 

0  <  tan'1  (-cot  a  sec  \p),  for  0  ^  ^  <  90°; 

9  >  tan'1  (-cot  a  sec  \p),  for  90°  <  ip  ^  180°. 

5.  RESULTS  AND  DISCUSSION  FOR  NON-NORMAL  DIRECTIONS  OF  E. 

Figures  10-17  show  escaping  secondary-electron  current  densities 
i  =  i(ar,^/,0,e)  when  -eE  is  not  normal  to  the  spacecraft  surface  (a  4  0°). 

Details  of  the  computations  of  these  results  are  the  same  as  those  given  in 
Secs.  2  and  3,  except  that  64  x  64  x  32  orbits  were  used  for  calculating  each 
value  of  i,  and  the  results  are  accurate  to  about  0.4%  or  better. 

In  Figs.  10  and  1 1,  t p  =  0°,  so  -eE,  B,  and  the  surface  normal  are  coplanar. 
The  E  x  B  drift  direction  is  therefore  tanpential  to  the  surface.  !r.  Fip.  10. 

O  o 

no  electrons  escape  (i  =  0)  when  -90°  ^  0  <  -60°  because  the  acceleration  of 


(4.3) 


(4.4) 


all  electrons  along  B  is  toward  the  surface,  as  discussed  in  Sec.  4.  For 
0  >  -60°,  this  acceleration  is  away  from  the  surface,  but  it  is  larger  for 
0  >  0°  than  for  0  <  0°,  so  increasing  the  electric  field  magnitude  e  increases 
electron  escape  more  for  0  >  0°.  For  0  >  -60°,  i  =  cos  0  when  e  =  0,  as  was 
the  case  for  -eE  normal  to  the  surface  (Fig.  5).  Some  care  is  needed  in 
defining  what  is  meant  by  the  case  "e  =  0".  For  -90°  ^  0  <  -60°,  i  =  0  in  the 
limit  e  -*■  0+  ,  but  in  the  limit  e  -+■  Q~,  the  electric  field  is  reversed,  so 
i  =  cos  0  in  this  limit,  and  i  therefore  has  a  discontinuity  at  6  =  0.  The 
discontinuity  is  reversed  for  0  >  -60°.  However,  we  are  interested  here 
primarily  in  cases  when  e  >  0  (negatively-charged  spacecraft  surfaces),  so  in 
this  work,  we  take  "e  =  0"  to  mean  the  limiting  case  e  ->  0-k  In  Fig.  1 1,  a 
has  oeen  increased  from  30°  to  60°,  so  the  effects  just  discussed  in  connec¬ 
tion  with  Fig.  10  are  seen  again,  but  more  strongly.  This  time,  escape  is 
suppressed  completely  for  -90°  ^  &  <  -30°. 

In  Fig.  12,  the  acceleration  of  electrons  along  B  is  toward  the  surface,  anc 
therefore  i  =  0,  for  0  <  tan_t(-V6)  =  -67.79°,  as  given  by  Eq.  (4.4).  Also,  we 
now  have  \p  4-  0°,  so  the  electric  force  vector  is  no  longer  in  the  same  Diane  as 
B  and  the  surface  normal.  As  a  result,  the  E  x  B  drift  now  has  a  nonzero 
normal  component.  This  decreases  electron  escape  for  6  <  0",  and  increases 
it  for  0  >  0°.  It  also  causes  the  escape  to  remain  nonzero  at  0  =  90°.  As 
before,  the  larger  outward  acceleration  along  B  also  increases  escape  for  0  > 
0°.  For  0  just  larger  than  67.79°,  we  see  that  escape  is  suppressed  almost 
completely  for  iarger  values  of  e;  this  is  because  the  inward  direction  of  the 
E  x  B  normal  component  causes  most  electrons  to  reimpact  the  surface 


during  the  first  gyroperiod  after  emission.  In  Fig.  13,  a  has  been  increased 
from  30°  to  60°,  with  consequent  enhancement  of  the  effects  just  discussed. 
We  new  have  i=0  for  0  <  tan_1(-V273)  =  -39.23°. 

In  rigs.  14  and  15,  \p  =  90°,  and  the  projection  of  -eE  along  B  is  away 
from  the  surface  for  all  0,  so  suppression  of  electron  escape  by  deceleration 
along  B  does  not  occur.  However,  for  0  close  to  -90°,  the  effect  of  the  inward 
direction  of  the  E  x  B  normal  component  overcomes  the  effect  of  the  outward 
acceleration  along  B,  especially  because  B  is  now  nearly  tangential  to  the 
surrace,  and  therefore  electron  escape  is  effectively  suppressed  for  larger 
values  of  e.  In  Fig.  15,  suppression  for  e  =  20  is  essentially  complete  over 
a  range  of  0  values  extending  more  than  30°  on  either  side  of  0  =  90°. 

in  Figs.  16  and  17,  we  have  \\t  =  135°,  and  the  effects  of  electron  deceler¬ 
ation  along  B  and  of  E  x  B  drift  now  suppress  electron  escape  at  oDposite  ends 
of  the  range  of  0.  From  Eq.  (4.4),  we  now  have  i  =  0  for  0  >  tan~'(V6)  = 
67.79°  and  0  >  tan''(V2/3)  =  39.23°,  respectively.  In  Fig.  16,  the  E  x  B 
effect  significantly  enhances  electron  escape  as  0  increases,  just  before  the 
deceleration  effect  cuts  it  off. 

The  results  shown  in  Fig.  5  and  Figs.  10-17  encompass,  albeit  rather 
sparsely,  the  entire  range  or  possible  directions  of  B  and  E  for  a  up  to  60°. 

To  see  this,  we  first  note  that  in  the  lmoortant  case  where  -eE  is  normal 


3f> 


to  the  surface,  the  i  values  for  -90°  ^0^0°  can  be  generated  from  those 
for  the  range  0°  ^  0  ^  90°,  which  is  covered  in  Fig.  5;  this  can  be  seen  by 
rotating  the  B  vector  in  Fig.  7  by  180°  about  the  z  axis.  Secondly,  the  i 
values  for  a  =  30°  and  60°  and  \p  =  180°,  225°,  270°,  and  315°  can  be 
generated  from  those  in  Figs.  10-17  by  rotating  the  B  and  -eE  vectors  in 
Fig.  7  together  about  the  z  axis  by  180°.  The  effect  of  this  is  to  increase  all 
the  i b  values  by  180°  and  also  to  reverse  the  sign  of  0.  Finally,  we  can  obtain 
the  i  values  for  cases  where  B  is  reversed  by  first  noting  that  reversal  of  B 
implies  reversal  of  both  the  £  coordinate  in  Fig.  7  (in  order  that  £  remain 
parallel  to  B)  and  the  q  coordinate  [in  order  that  the  (^y,q)  axes  remain 

right-handed].  The  quantities  sin  0,  cos  0,  u^,  uq^,  f,  e^,  and  in  Eq.  (2.8) 
will  then  all  reverse.  To  keep  z  =  z(r)  in  Eq.  (2.8)  unchanged,  we  require  also 

that  n  reverse,  and  we  therefore  require  that  e  and  u  also  be  reversed. 

^  y  oy 

However,  the  emitted  velocity  distribution  is  symmetric  in  u  ,  and  reversing 
6^  involves  replacing  the  angle  \\i  by  360°  -  \p,  and  this  replacement  gives  back 
the  same  set  of  \p  values  for  which  our  computations  already  give  i. 


We  can  summarize  the  results  in  Figs.  10-17  by  noting  that  when  a  ^  0, 
two  new  mechanisms,  which  were  not  present  when  -eE  was  normal  to  the 
surface  (Fig.  5),  car.  suppress  electron  escape.  These  are:  an  inward  normal 
component  of  E  x  B,  if  e  is  strong  enough,  and  a  decelerating  projection  of 
-eE  along  B,  for  any  e  >  0.  These  may  act  at  the  same  end  or  at  opposite  ends 
of  the  range  of  magnetic-field  directions  -90°^  0  ^  90°.  These  mechanisms 
can  greatly  enlarge  the  range  of  surface  orientations  for  which  escape  is 
suppressed. 


6.  CALCULATION  OF  SECONDARY-ELECTRON  DENSITIES 


Once  the  secondary-electron  escape  fluxes  are  known  (Secs.  3-5),  a  simple, 

inexpensive,  approximate  calculation  of  their  space-charge  density  distribution 

can  be  set  up.  The  proposed  method  is  as  follows:  (1)  ignore  the  gyromotion 

of  the  secondary  electrons  once  they  have  escaped.  Their  motion  then  involves: 

(a)  an  acceleration  along  magnetic  field  lines,  of  amount  -(e/m)E»B/B  (b)  a 

drift  motion  of  velocity  E  x  B/B2  across  magnetic  fieid  lines.  (2)  Integrate 

enough  of  the  trajectories  defined  by  this  motion  (i.e.  their  guiding-center 

trajectories)  to  define  trajectory  tubes  whose  cross-section  at  any  point  can  be 

calculated  with  sufficient  accuracy;  the  method  described  by  Laframboise  et  al 

(1982b,  Sec.  7),  can  be  used  to  calculate  the  area  of  a  trajectory  tube  without 

reference  to  neighbouring  trajectories.  (3)  Calculate  their  space-charge 

density  n(r  )  at  any  point  by  (a)  ignoring  the  "thermal"  spread  of  their 

velocities  (b)  then  invoking  the  fact  that  their  density  x  their  velocity  [as  given 

by  the  orbit  integration  mentioned  in  (2)),  x  the  cross-sectional  area  A(r  )  of 

the  trajectory  tube  (which  must  be  calculated  in  a  plane  J.  the  trajectory)  at 

the  point  r  in  question,  =  a  constant  (whose  value  is  given  by  the  initial 

conditions  at  the  point  on  the  spacecraft  where  the  trajectory  originates)  (c) 
finding  their  velocity  at  the  point  in  question  by  using  energy  conservation, 

together  with  the  values  of  electric  potential  <p(r )  and  <pQ  at  that  point  and  the 

emission  point,  and  their  assumed  velocity  vq  at  the  emission  point.  The 

result  is: 

r'r'  =  °  v  A  /  [A(r)vv  2  +  (2e/mi  [<£(r  '  -  &>  ]  1 

o  o  o  (  -  O  vOi  J 


(6.1) 


where  nQvo  is  the  escaping  flux  calculated  in  Secs.  3-5.  At  most  positions, 
n(r)  will  be  insensitive  to  the  precise  value  assumed  for  vq2;  assuming  that 

vq  =  the  one-sided  thermal  speed  (2kT/7rm)^  will  suffice  for  most  purposes. 
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TABLE  1 

Values  of  the  ratio  i  =  I/Iq  of  escaping  to  emitted  flux,  for  various  values 
of  0,  the  angle  (in  degrees)  between  the  surface  normal  and  the  magnetic  field 
direction,  and  e,  the  nondimensional  repelling  electric  field  strength.  These 
two  quantities  appear  in  the  table  as  THETA  and  EPS,  respectively.  These 
results  are  accurate  to  within  about  0.2%  or  better.  The  electric  field  is 


normal  to  the  surface. 


FIGURE  CAPTIONS 


Figure  1.  Effect  of  surface  orientation  on  escape  of  emitted  electrons.  In  (a), 
the  spacecraft  surface  is  perpendicular  to  the  magnetic  field  B  ,  and  the 
emitted  electrons,  which  experience  an  electric  force  -eE  directed  away 
from  the  surface,  all  escape.  In  (b),  the  spacecraft  surface  is  nearly 
parallel  to  B,  and  almost  all  of  the  emitted  electrons  return  to  the  surface, 
even  though  they  still  experience  an  electric  force  directed  away  from  it. 
Note  that  the  component  of  E  perpendicular  to  B  results  only  in  an  E  x  B 
drift  parallel  to  the  surface. 

Figure  2.  Spacecraft  simultaneously  in  a  collisionless  ion  flow  and  a  magnetic 
field  B. 

Figure  3.  Coordinate  system  for  calculating  electron  escape  fluxes  when  E  is 
perpendicular  to  the  spacecraft  surface.  The  y-coordinate  (not  shown)  is 
directed  into  the  plane  of  the  Figure. 

Figure  4.  Example  of  an  electron  orbit  having  zero  initial  velocity.  The 

magnetic  field  B  is  parallel  to  the  (x,z)  plane,  and  makes  an  angie  Q  =  75° 
with  the  z  axis,  e  =  1 .  Three  gyroperiods  of  the  orbit  (0  ^  r  ^  6rr)  are 
shown. 
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Figure  5.  Ratio  i  =  I/I  of  escaping  to  emitted  secondary  electron  flux,  as  a 
function  of  the  angle  0  between  the  surface  normal  and  the  magnetic  field 
direction,  for  various  values  of  the  repelling  electric  field  strength  param¬ 
eter  e  =  (E/'B)  (7rm/'2kT)*.  The  eiectnc  field  is  normal  to  the  surface. 

The  result  for  e  =  0  is  given  by  l  =  cos  0.  Realistic  values  of  e  for  Shuttle 
high-voltage  charging  conditions  are  in  the  range  30  $  t  u  120  (Sec.  2)! 

Figure  6.  Electron  orbit  for  e  =  0,  fictitiously  extended  so  as  to  pass 
through  the  surface  and  re-emerge  from  it. 

Figure  7.  Problem  geometry  when  potential  also  varies  along  surface. 

Figure  8.  (a)  Typical  orbit  of  emitted  electron  when  the  electric  force 

-eE  on  it  has  an  outward  normal  component,  but  the  pro  lection  of  -eE 
along  B  has  an  inward  normal  component  (b)  analogous  situation 
involving  sailboat  tacking  into  wind. 

Figure  9.  Dependence  of  secondary-electron  escape  ana  surface  currents  on 
electric  field  direction  at  surface. 

Figure  10.  Ratio  i  =  I/I  of  escaping  to  emitted  secondary  electron  flux,  as  a 
function  of  the  angle  0  between  the  surface  normal  and  the  magnetic-field 
direction,  for  various  values  of  the  electric  field  strength  parameter 
e  =  (E/B)  (7rm/2kT)  Same  as  Figure  5,  except  that  the  electric  force 
vector  -eE  IS  no  longer  normal  to  the  surface  (a  is  nonzero). 

Figure  1  1.  Same  as  figure  10,  except  that  -eE  is  tilted  further  away  from  the 
surface  normal  (a  =  60°). 


Figure  12.  Same  as  Figure  10,  except  that  -eE  is  no  longer  in  the  same  plane 
as  the  surface  normal  and  the  magnetic  field  vector  (\p  is  nonzero). 

Figure  13.  Same  as  Figure  12,  except  that  -eE  is  tilted  further  away  from  the 
surface  normal. 

Figure  14.  Same  as  Figures  10  and  12,  except  that  i//  =  90°. 


Figure  15.  Same  as  Figure  14,  except  that  -eE  is  tilted  further  away  from  the 
surface  normal. 


figure  16.  Same  as  Figures  10,  12,  and  14,  except  that  =  135°. 

Figure  17.  Same  as  Figure  16,  except  that  -eE  is  tilled  further  away  from  the 
surface  normal. 
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Figure  1.  Effect  of  surface  orientation  on  escape  of  emitted  electrons.  In 
(a),  the  spacecraft  surface  is  perpendicular  to  the  magnetic  field  B,  and  the 
emitted  electrons,  which  experience  an  electric  force  -eE  directed  away  from 
the  surface,  ail  escaDe.  In  (b),  the  spacecraft  surface  is  nearly  parallel  to  B, 
and  almost  all  of  the  emitted  electrons  return  to  the  surface,  even  though  they 
still  experience  an  electric  force  directed  away  from  it.  Note  that  the 
comDorent  of  E  perpendicular  to  B  results  only  in  an  E  x  B  drift  parallel  to 
the  surface. 
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2.  Spacecraft  simultaneously  in  a  collisionless  ion  flow  and  a  magnetic 
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Figure  3.  Coordinate  system  for  calculating  electron  escape  fluxes  when  E  is 
perpendicular  to  the  spacecraft  surface.  The  y-coordinate  ;rct  snown)  is 
directed  into  the  plane  of  the  Figure. 
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gure  A.  Exarr.p.e  of  an  electron  orbit  having  zero  initial  velocity.  The 
magnetic  fieid  B  is  parallel  to  the  (x,z)  plane,  and  makes  an  angle  0  =  7, 
jvith  the  z  axis.  6  =  1.  Three  gyroperiods  of  the  orbit  (0  ^  r  ^  6 7 r)  are 
shown. 
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ure  5.  Ratio  ;  =  !/I0  of  escaping  to  emitted  secondary-electron  flux,  as  a 
function  of  the  angie  0  between  the  surface  normal  and  the  magnetic  field 
direction,  for  various  values  of  the  repelling  electric  field  strength  param¬ 
eter  e  -  (E/B)  (Trm/2kT)J  .  The  electric  field  is  normal  to  the  surface. 

The  result  for  r  =  0  is  given  by  i  =  cos  0.  Realistic  values  of  e  for  Shuttle 
high-voltage  charging  conditions  are  in  the  range  30  $  f  $  120  (Sec.  2)! 
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Figure  7.  Problem  geometry  when  potential  also  varies  along  surface. 
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Figure  8.  (a)  Typical  orbit  of  emitted  electron  when  the  electric  force 
-eE  on  it  has  an  outward  normal  component,  but  the  projection  of  -eE 
along  B  has  an  inward  normal  component  (b)  analogous  situation  involv¬ 
ing  sailboat  tacking  into  wind. 
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igure  9.  Dependence  of  secondary-electron  escape  and  surface  currents  on 
electric  field  direction  at  surface. 
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Figure  10.  Ratio  i  =  I/I  of  escaping  to  emitted  secondary-electron  flux, 
as  a  function  of  the  angle  9  between  the  surface  normal  and  the 
magnetic-field  direction,  for  various  values  of  the  electric  field  strength 
parameter  e  =  (E/B)  (7rm/2kT) *.  Same  as  Figure  5,  except  that  the  elec¬ 
tric  force  vector  -eE  is  no  longer  normal  to  the  surface  (a  is  nonzero). 
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igure  17.  Same  as  Figure  16,  except  that  -eE  is  tilted  further  away  from  the 
surface  normal. 
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